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Abstract. Each spring a layer of small particles forms between 20 and 30 km in the polar regions. Results are presented 
from a 2D microphysical model of sulfate aerosol, which provide the first self-consistent explanation of the observed ”CN 
layer.” Photochemical conversion of sulfuric acid to SO, in the upper stratosphere and mesosphere is necessary for this layer 
to form. Recent laboratory measurements of H 2 S0 4 and SO, photolysis rates are consistent with such conversion, though an 
additional source of SO, may be required. Nucleation throughout the polar winter extends the top of the aerosol layer to 
higher altitudes, despite strong downward transport of ambient air. This finding may be important to heterogeneous 
chemistry at the top of the aerosol layer in polar winter and spring. 
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Introduction 

Rosen and Hofmann [1983] first noted an anomalous enhancement in springtime particle concentration near 30 km 
altitude at mid-latitudes in air of recent Arctic origin. Similar enhancements were observed from September to November 
over Antarctica [ Hofmann and Rosen, 1985; Hofmann, 1988; Hofmann et al., 1988; Hofmann et al., 1989; Wilson et ai, 
1989; Hofmann, 1990; Hofmann and Deshler, 1991], and in January and February over the Arctic [Hofmann, 1990; Wilson et 
al., 1990]. Because the concentration is dominated by the smallest particles, known as “condensation nuclei,” the polar 
anomaly has become known as a “CN layer.” 

Antarctic in situ measurements indicate modest particle nucleation at the top of the sulfate layer before the arrival of 
sunlight, with dramatic increases producing a well-defined layer when sunlight arrives. Balloon-borne measurements from 
McMurdo Station in late August 1987 and 1988 show particle concentrations increasing to a maximum of -20 particles/cm' 
near 25 km [Hofmann et ai, 1989; Hofmann, 1990]. Because particle concentrations should otherwise decrease with air 
density (as is observed in polar summer), these observations suggest that nucleation occurs between 20 and 30 km in polar 
night. Observations made in each austral spring from 1983 to 1989 are consistent with a well-defined layer of 100-150 cm ' 
forming each September between 22 and 28 km, with the altitude of the peak concentration descending through October to 
near 20-22 km in early November [Hofmann and Rosen, 1985; Hofmann, 1988; Hofmann et ai, 1988; Hofmann et al., 1989; 
Hofmann, 1990; Hofmann and Deshler, 1991], In 1993, the year following the eruption of Mt. Pinatubo, particle 
concentrations reached 50 cm 5 by late August, and peaked at about 200 cm' 3 in early September [Deshler et ai, 1994]. 

The CN layer was first discovered in Northern Hemisphere measurements, and measurements clearly show well- 
defined CN events in the arctic, albeit with less consistency than the Antarctic. Observations over Kiruna. Sweden (68°N) 



show a layer of 40/cm 3 between 22.5 and 26 km on January 30, 1989 [ Hofmann , 1990]. But four profiles over Kiruna in 
January and February of 1992 show no definite layer, with CN concentrations constant or increasing slightly with altitude 
[Deshler, 1994]. The increasing particle mixing ratio with altitude in these profiles indicates that some nucleation occurred at 
the top of the layer, comparable to the weak nucleation indicated in Antarctic winter before sunrise. The variability in these 
arctic measurements are likely a product of the greater latitudinal excursions of vortex air in the Arctic as opposed to the 
Antarctic, as well as the lack of measurements well above the Arctic Circle. Trajectory analyses of the air sampled over 
Kiruna in 1992 indicated that it had recently been as far south as 50°N, and may have mixed with mid-latitude air. Such 
excursions and mixing explain the observation of the arctic CN layer over Laramie, Wyoming (41°N) in late winter and early 
spring of 1982, 1983, and 1986 [Rosen and Hofmann , 1983; Hofmann el ai , 1985; Hofmann et ai y 1988]. Because the 
Antarctic vortex is more circumpolar, leading to less variation in latitude and sunlit history of air parcels, and because 
Antarctic observations of the CN layer are more extensive and consistent than those in the Arctic, the focus of this paper is 
the Antarctic CN layer. It should be noted, however, that while the dynamical variability makes modeling the timing of the 
CN layer more complicated in the Arctic, the basic theory we present applies there as well. 

Theories 

Rosen and Hofmann [1983] suggested two possible origins for the CN layer. First, sulfur-bearing gases might 
maximize in winter, due to lack of sunlight available to oxidize them to sulfuric acid. These gases might dissociate rapidly 
with the return of sunlight, OH, and O in spring. Their second theory relies on very rapid drops in temperature observed in 
the Arctic to create very large supersaturations. In the absence of large, existing particles on which to condense, the excess of 
ambient vapor would nucleate many new, small particles. 

Building on the nucleation rate studies of Hamill et ai [1990], Zhao et ai [1995] investigated the CN layer with a 1- 
dimensional microphysical model. They questioned whether sufficient sulfuric acid is present in polar winter, given the lack 
of photochemical production there, to initiate the temperature-dependent mechanism of Rosen and Hofmann [1983]. They 
also calculated that OCS oxidation can not provide sufficient sulfuric acid vapor to allow new particles to grow to observable 
size, as suggested by Oppenheimer [1987]. They postulated that the “extra” sulfur required for the formation of the layer can 
be provided by SO, descending from the mesosphere. The source of this SO, was hypothesized to be H 2 S0 4 photolysis at 
lower latitudes. By setting upper boundary conditions of 100 pptv for both SO, and H,S0 4 at 60 km and allowing for vertical 
subsidence, Zhao et ai were able to simulate the formation of a CN layer similar to that which is observed. The CN 
concentration was found to be much more sensitive to the SO, concentration than the H,S0 4 concentration. SO, 
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concentrations increase during the polar night as advection carries SO : -rich air downward. During descent gas-phase sulfuric 
acid is lost on existing particles, and S0 2 remains as the dominant gaseous sulfur reservoir until the return of sunlight in 
spring. Zhao et al. did not model the production of SO, based on measured photolysis rates as we do here. 

Photochemistry of S0 2 

We present here calculations from a new microphysical model, which has been incorporated into the Garcia- 
Solomon two-dimensional dynamical/chemical model. The model spans 56 pressure levels from 2 to 112 km above sea 
level, and 36 latitude bins from 89.5°S to 89.5°N. Details of the photochemistry and dynamics of the model are described in 
Garcia et al. [1992] and Garcia and Solomon [1994]. We have added both sulfur oxidation chemistry and particle 
microphysics to this model. 

OCS, tropopause SO,, and tropopause aerosol are the sources of stratospheric aerosol in the model. Oxidation of 
gas-phase sulfur species is calculated according to the recommendations of DeMore et al. [1994] except as noted below. The 
mixing ratio for OCS is maintained at 510 pptv at the lower boundary of the model (750 mbar). OCS is the dominant source 
of sulfate at the top of the stratospheric aerosol layer. In the lower stratosphere, most of the sulfate comes from tropospheric 
aerosol and SO,, which are maintained at constant mixing ratios of 0.27 ppbm of SO/ and 80 pptv of SO, at the model's 
tropical tropopause. Tropospheric aerosol is represented by a bimodal lognormal distribution, with 30% of the mass in a 0.8- 
jim radius mode, and the majority of particles distributed about 0.015pm. The details of these tropospheric sources of sulfur 
do not affect the calculation of the polar CN layer. 

As temperatures increase with altitude in the stratosphere, sulfate aerosol evaporates, producing sulfuric acid vapor. 
This vapor is thought to photolyze to SO, at the short wavelengths of sunlight present in the upper stratosphere and 
mesosphere. Rinsland et al. [1995] concluded from observations of SO, in the upper stratosphere and mesosphere and from 
2-dimensional model calculations, that sulfuric acid vapor should photolyze at rates roughly 0.3 times that of HCl. We 
present this assumption as model case A (see table 1). Recently, Burkholder et al. [1998] measured no photolytic cross 
sections for H,SO, vapor greater than their upper limit of 10' 21 cm 2 , for wavelengths greater than 195 nm. Our model case B 
sets the cross sections to that upper limit of 10' 21 cm 2 in the 179-226 nm region where HCl absorbs, representing maximum 
photolysis rates consistent with the findings of Burkholder et al. Case C assumes that the cross section is 10 21 cm 2 at 195 nm, 
with the same wavelength dependence as HCl (decreasing with increasing wavelength). This assumption sets the cross 
sections to 0.01 1 times those of HCl. Case D assumes that H,S0 4 vapor does not photolyze at all. 
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Burkholder and McKeen [1997] recently measured S0 3 photolysis cross sections, which we have included in the 
model. This reaction is an important source of S0 2 near the top of the sulfate layer, where S0 3 photolyzes faster than it 
oxidizes to H„S0 4 . The oxidation of S0 3 proceeds through the production of an S0 3 *H 2 0 adduct, resulting in a pressure- 
dependent effective reaction rate [Lovejoy et at, 1996]. The timescale for oxidation increases exponentially with altitude, 
from a few seconds near the tropopause to a few days at 40 km. Photolysis becomes faster with altitude, overtaking 
oxidation as the dominant loss for S0 3 near 40 km. Model cases A through D include S0 3 photolysis. We test the 
importance of the new S0 3 photolysis reaction to the formation of the CN layer with case E and F, which are the same as 
cases B and D, respectively, but without S0 3 photolysis. 

Calculated SO, mixing ratios for cases B, D, E, and F are shown in Figure I, revealing the effects of photolysis of 
H,S0 4 , S0 3 , and SO,. The figures show SO, on June 1 building up to more than 170 pptv in case B. Mixing ratios are 
greatest in polar winter, where destruction by photolysis is shut off, and SO thus formed reverts to SO,. 

Particle Microphysics 

The microphysical model includes classical binary homogeneous nucleation; coagulation; condensation and 
evaporation of water, sulfuric acid, and nitric acid; transport with ambient air parcels; and sedimentation [Mills, 1996]. The 
model simulates the formation of Type lb (liquid ternary solutions) and Type II (water ice) polar stratospheric clouds, but 
does not deal with Type la (nitric acid hydrate) PSCs. Aerosol particles are divided into 45 size classes ranging from -0.5 
nm to -10 pm. 

Calculation of binary homogeneous nucleation rates is highly uncertain. Wyslouzil et ai [1991] present laboratory 
measurements of nucleation rates, showing discrepancies of ±10 orders of magnitude with classical theory. Weber et ai 
[1996] conclude from atmospheric measurements that particle number varies as the square of sulfuric acid concentrations, 
rather than the classically predicted exponential response. Brock et ai [1995] point out, however, that despite these 
uncertainties, classical nucleation theory correctly predicts the regions in the atmosphere in which large numbers of particles 
are observed. After experimenting with alternative nucleation parameterizations based on the findings of Wyslouzil et ai and 
Weber et ai , we conclude that classical nucleation theory provides the best correspondence with observed aerosol behavior in 
the stratosphere. 

Figure 2 shows the calculated concentrations of particles larger than 0.01 pm, the lower size limit of balloon-borne 
CN counters, from the 6 model cases described in Table 1. Cases A and B reproduce well the observed particle 
concentrations in both winter and spring. In these cases, particle production is evident throughout the winter, with dramatic 
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increases to 80-120 particles/cm 1 in September. The CN layer descends from 30 to 22 km as the austral spring progresses. 
Measurements made repeatedly in the spring of a given year show the CN layer to be descending in a similar manner 
[Hofmann et al y 1988; Hofmann et al , 1989; Deshler et al . , 1994]. Although the observed CN layer sometimes initiates as 
low as 26 km, these observations also indicate the presence of enhanced CN concentrations at 30 km. The calculations 
diverge from observations in the other model cases, when either H 2 S0 4 or S0 3 photolysis rates are decreased. 

Figure 3 shows the seasonal cycle of temperature, nucleation, and particle concentration for case B. There is a 
small, but significant rate of nucleation near 30 km as early as April which continues throughout the winter, as cold 
temperatures (below 220 K) reach higher altitudes. Prior to July, at any given altitude nucleation returns to negligible rates 
within several weeks of its initiation, despite continued decreases in temperature, as aerosol formation and growth deplete 
ambient sulfuric acid vapor. The winter phase of nucleation produces small concentrations (10-50 /cm 3 ) ot particles, most of 
which grow to observable size at the 30 km maximum. Observations at McMurdo Station (78°S) show particle 
concentrations increasing slightly from about 10 /cm 3 at 16 km to about 20/cm 3 at 26 km in late August, before the well- 
defined CN layer forms [ Hofmann et al , 1989; Hofmann , 1990]. Our calculations suggest similar concentrations of 
observable particles at 26 km. Aerosol is not suppressed by the strong descent prevalent in polar winter, but instead is 
present at higher altitudes than in summer due to nucleation. If the only source of particles were at lower altitudes (i.e. the 
tropopause), subsidence would push down the top of the layer in polar stratospheric winter. 

About a month after temperatures reach their minimum below 190K in late June and early July, particle nucleation 
increases dramatically at all altitudes below about 36 km. Particle concentrations respond rapidly, peaking at over 107cm 3 
between 22 and 30 km. This event coincides with the return of sunlight in late July, and continues until temperatures 
increase to greater than 200 K at the beginning of October. A layer of observable-sized particles forms above 30 km in early 
September, descending in altitude and increasing in concentration until a maximum of about 80/cm 3 is found near 26 km in 
late October. This maximum in observable particles occurs nearly 2 months after the peak in total particle concentration. 

The delay between the minimum in temperature and the maximum in calculated nucleation, and the subsequent lag 
between concentrations of total particles and observable particles depend critically on the descent of S0 2 from above. Figure 
4 shows the annual cycles of various calculated gas-phase species and aerosol parameters at 26 km and 74°S for cases B, D, 
and E. The arrival of polar night at 74°S is evident from the dramatic drop in OH as photolytic production shuts off at the 
end of May. Air descends from the mesosphere to the stratosphere throughout polar winter. In case B, the air carries with it 
mixing ratios of SO : as high as 160 pptv at some altitudes (see Figure 1). Liquid sulfate mass, shown in Figure 4, decreases 
by sedimentation and descent throughout the austral winter in all three cases. 
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In mid-May, the disappearance of OH cuts off all production of H 2 S0 4 vapor, while the rate of loss due to 
condensation increases slightly on newly nucleated particles descending from above 30 km (see Figure 3). Loss of 
production and enhancement of the sink causes H 2 S0 4 mixing ratios to plummet by 2-4 orders of magnitude in cases B, D, 
and E. The effect is greatest in case D because more H,S0 4 vapor is available for nucleation above 30 km in May, resulting 
in greater surface areas and faster condensation. 

Sulfuric acid vapor recovers in July, coincident with the return of sunlight and OH. At that time production of 
H,S0 4 vapor from SO, resumes (quite rapidly in cases B and E) as SO, has increased to more than 50 times its summertime 
minimum. The increase in sulfuric acid vapor concentration by two orders of magnitude in July, while cold temperatures 
continue to suppress equilibrium vapor pressures, generates significant nucleation rates, especially in cases B and E. 
Sustained over the months of August and September, the nucleation event produces particle concentrations in excess of 10 4 
/cm 3 in these cases (see Figure 3). 

By late September these particles at 26 km grow to observable size in case B. The growth in particle size coincides 
with a dramatic increase in liquid sulfate mass, indicating that this growth is not due merely to coagulation. The steady 
increase in OH with increasing sunlight as spring progresses hastens the conversion of SO, to sulfuric acid, which condenses 
on the recently nucleated particles, rather than nucleating new ones. Nucleation shuts off in early October in favor of 
condensation due to rising temperatures. At 26 km in polar winter, homogeneous nucleation rates calculated from classical 
theory decrease by about 6 orders of magnitude as temperature increases from 180K in late September to 200K in early 
October. Condensation rates have a much lesser dependence on temperature. 

Comparison of cases B, D, and E in Figure 4 reveals the importance of both SO, and H 2 S0 4 photolysis to the 
production of a CN layer resembling that observed. H 2 S0 4 photolysis provides an essential mechanism for suppressing 
homogeneous nucleation in winter and conserving gas-phase sulfur as SO,, which becomes available again as H,S0 4 in the 
spring to fuel nucleation. Although Burkholder et al.'s measurements reveal no evidence for H,S0 4 photolysis, their upper 
limit produces enough SO,, in conjunction with SO, photolysis, to create a CN layer similar to that observed. SO, photolysis 
prevents SO, from oxidizing back to H,S0 4 . At 40 km SO, photolysis can increase descending SO, mixing ratios by 50 pptv 
(see Figure 1), and at 26 km it increases the peak SO, mixing ratio seen in September by a factor of 2.25 (see Figure 4). 

Figure 5 shows the effect of the CN layer on aerosol surface areas at 74°S in each of the six model cases. In all 
cases, surface areas between 30 and 35 km are enhanced in the winter and spring relative to summer and autumn. H,S0 4 
photolysis in cases A and B results in springtime surface areas enhanced by as much as a factor of 2 over cases D and F, in 
which H,S0 4 photolysis is turned off. Polar winter nucleation and the springtime CN layer may affect heterogeneous 
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reactions between sulfate aerosol and various gas-phase species on a seasonal basis. In all cases, surface areas of 0.1 p.m7cm' 
or greater extend up to 36 km throughout winter and spring, and to less than 30 km in summer and autumn. Depending on 
reaction probabilities and ambient concentrations of gas-phase species, such surface areas could perturb gas-phase species via 
heterogeneous chemistry. In the absence of new particle nucleation above 30 km, dynamical descent in the polar winter 
would suppress the sulfate aerosol layer, reducing rather than increasing the altitude of the top of the layer relative to 
summer. 

Conclusions 

The 2D microphysical model of sulfate aerosol that we have developed provides the first self-consistent explanation 
of the CN layer observed in Antarctic spring. Although dynamics complicate the Arctic situation, the same theory applies to 
enhanced CN observed there [ Hofmann , 1990; Wilson et ai, 1990; Deshler , 1994], 

Binary homogeneous nucleation and condensation following rapid production of sulfuric acid vapor from SO, in 
spring is responsible for the appearance of a CN layer, as suggested by Zhao et al. [1995]. Photolysis of S0 3 and H,S0 4 are 
proposed as the mechanism for producing S0 2 in the upper stratosphere and mesosphere; downward transport in polar winter 
and spring then bring the SO, into the lower stratosphere. The model demonstrates that photolysis of S0 3 and H,S0 4 
calculated within the constraints of recent laboratory measurements, can partition sulfur as SO, in the upper stratosphere and 
mesosphere to produce a springtime CN layer much like that observed. However, when a wavelength dependence in the 
H,S0 4 cross sections matching that of HC1 further constrains photolysis, the calculated CN layer has substantially smaller 
concentrations than that observed. An additional source for S0 2 in the upper stratosphere and mesosphere may be indicated 
to explain the CN layer, as well as observations of SO, at these altitudes. 

The model also suggests that new particles nucleate and grow throughout the winter above 30 km, at a slower rate 
than that which produces springtime CN layer. Thus aerosol is present at higher altitudes in polar winter than in summer, 
despite the subsidence of polar air which would otherwise push the top of the aerosol layer downward. Concentrations of 
new particles thus formed are consistent with in situ measurements made in polar winter. These particles may be important 
to heterogeneous chemistry at the top of the aerosol layer in polar winter and spring. 

Aknowledgments. This work was supported in part by the NASA SAGE II Science Team through grant L68687D 
and NASA's Atmospheric Chemistry Program NAG 1-2021. 
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Tables 


Case 

H 2 S0 4 cross sections * 

A 

^H2S04 — 3 *^HC: 

B 


C 


D 

W’ i-nn ■ .- 

Vn2SO< u,v/ ' 

E 

^H2S04~“ 1 P. ^ 

F 

^H2S04 0.0 


: . ^'J^fS'SOiphotolysis 

Rinsland fit to observations On 


„=10' JI cm 2 . Burkholder et al. upper limit ! V; ' On 

Upper limit with HC1 tail On 

Burkholder et al. upper limit Off 

No H,S0 4 photolysis Off 


Table 1. Model cases, varying H 2 S0 4 and SO, photolysis. 


Figure Captions 

Figure 1. SO, mixing ratios (pptv) for four model cases representing various assumptions about H,S0 4 and SO, photolysis. 
Calculations are shown for noon, Junel. 

Figure 2. Calculated concentrations (/cm 5 ) of particles large enough to be observed (r>0.01|im) in six model cases, varying 
photolysis rates for H 2 S0 4 and SO,. Results are presented as annual cycles at 74°S, averaged over 5-day periods for each 
altitude gridpoint (2 km resolution). In cases A and B, CN layers are calculated which resemble those observed. SO, 
photolysis is included in cases A through D. Comparison of cases B and E reveals the importance of SO, photolysis to the 
concentration of particles in the CN layer. 

Figure 3. Model calculations of (a) temperatures, (b) homogeneous binary nucleation rates, (c) total particle concentrations, 
and (d) particles of radius greater than 0.01pm at 74°S in model case B. A delay of about one month is evident between the 
periods of minimum temperature and maximum nucleation. A similar delay occurs between the maximum in total particle 
concentration and that in particles large enough to be observed. 

Figure 4. Variation in OH, S02, sulfate mass, aerosol surface area, sulfuric acid vapor, sulfuric acid equilibrium vapor 
pressure, nucleation rate, and particles greater than 0.01 micron radius at 26 1cm, 74°S for model cases B, D, and E (see table 
l and text). 

Figure 5* Calculated CN layer surface areas (pmVcm 3 ) in six model cases at 74°S. 
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